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SUMMARY 


The  CRM  differential  equations  are  stated  in  simplified 
form,  and  refined  values  for  model  parameters  and  initial  condi> 
lions  are  presented  and  discussed. 

Detailed  comparisons  of  simulation  results  with  observed 
data  are  presented  over  the  yield  range  from  0.5  T to  15  MT,  with 
emphasis  on  the  range  10  T to  15  MT.  Comparisons  are  mostly  for 
stabilized  cloud  dimensions,  but  some  comparisons  of  temporal 
cloud  development  are  included.  Considering  the  scatter  and  un- 
certainty of  observed  data,  the  comparisons  indicate  excellent 
agreement  for  the  most  part,  though  there  may  be  a tendency 
to  overpredict  base  heights  for  very  low  yield  shots  (see  Figs.  1-6). 

Compared  with  a simple  power  function  of  yield  fitted  to 
the  observed  data  by  least  squares  (eq.  (13)),  CRM  results  yield  a 
59*  reduction  in  the  variance  of  1o9(2-j-)^g1j.  - ^09(zj)Q|jg* 

Non-linearity  of  the  model  makes  it  very  difficult  to 
estimate  a priori  the  effects  of  changes  in:  model  structure,  para- 
meter values,  and  initial  conditions.  Results  are  sensitive  to 
changes  in  cloud  shape,  the  turbulence  parameter  k^,  and  the  en- 
trainment parameter  u.  Furthermore,  these  interact  in  complex 
ways.  One  must  proceed  with  caution  in  making  changes,  and  effects 
of  changes  must  be  assessed  via  thorough  study  of  simulation  re-  • 
suits  across  the  complete  range  of  yields. 

Fundamental  deficiencies  in  the  model,  which  are  es- 
pecially acute  at  extremely  high  yields,  are  pointed  out.  Extra- 
polation of  usage  to  the  yield  range  20-100MT  results  in  apparent 
anomalous  results,  which  are  explained  as  resulting  fn.m  the  de- 
ficiencies. Use  of  the  model  for  such  high  yields  is  question- 
able, especially  since  observations  are  not  available  for  model 
calibration  and  validation. 
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1.  INTRODUCTION 

Prediction  of  local  fallout  transport  and  deposition  Is 
critically  dependent  on  the  height  and  dimensions  of  the  stabilized 
nuclear  cloud.  Most  prediction  calculations  begin  with  stabilized 
clouds  whose  dimensions  are  similar  to  those  produced  by  test  shots 
at  the  Nevada  Test  Site  and  the  Marshall  Islands.  This  provides  no 
flexibility  to  account  for  effects  of  atmospheres  that  have  dif- 
ferent stability  structures  from  those  typical  of  the  test  locations. 
The  OELFIC  Cloud  Rise  Module  (CRM)  uses  a dynamic  model  of  cloud 
rise  that  adequately  accounts  for  atmospheric  stability,  as  well 
as  for  many  other  variables  that  may  be  significant. 

The  CRM  Is  based  on  the  water-surface-burst  cloud  model 

(1  21 

developed  at  the  Naval  Radiological  Defense  Laboratory  (NRDL)^  ' . 

Details  of  the  mathematical  basis  of  the  model  and  of  the  CRM  code 

f3l 

are  given  by  Norment  and  Woolf' 

Various  Individuals  have  worked  on  application  of  the 
model  to  land-fallout  prediction  over  the  span  of  ten  years  since 
it  was  first  adapted  for  use  In  DELFIC.  However,  the  work  reported 
here  Is  the  first  really  comprehensive,  thorough  calibration  and 
validation  study  that  has  been  done. 

The  CRM  model  is  refined  and  its  strengths  and  weak- 
nesses are  explored  and  described. 
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2.  THE  DELFIC  CLOUO  RISE  MODULE 


2.1  DESCRIPTION 

The  DELFIC  CRM  is  a dynamic,  one-dimensional,  entrain- 
ing buoble  model  of  nuclear  cloud  rise.  It  consists  of  a set  of 
coupled  ordinary  differential  equations  that  represent  conserva- 
tion of  momentum,  mass,  heat  and  turbulent  kinetic  energy.  The 
nuclear  cloud  is  defined  in  terms  of:  vertical  coordinate  of  its 
center  (the  cloud  is  in  some  respects  treated  as  a point)  , cloud 
volume,  average  temperature,  average  turbulent  energy  density, 
and  the  masses  of  its  constituents:  air,  soil,  weapon  debris, 
water  vapor  and  condensed  water.  Cloud  properties  and  contents 
are  taken  to  be  uniform  over  the  cloud  volume. 

Initial  conditions  are  specified  at  approximately  the 
time  the  fireball  reaches  pressure  equilibrium  with  the  atmosphere. 
Atmospheric  conditions  (vertical  profiles  of  pressure,  temperature 
and  relative  humidity)  are  accepted  by  the  CRM  in  tabular  form. 

Effects  of  initial  conditions  and  atmospheric  stability 
on  the  cloud  rise  are  accounted  for  because  of  the  dynamic  nature 
of  the  model.  Atmospheric  stability  is  especially  important  in 
determining  the  stabilized  height  of  the  cloud,  which,  in  turn,  is 
critical  in  determining  the  distribution  of  local  fallout  on  the 
ground. 

Thr  differential  equations  are  solved  by  a fourth-order 
Runge-Kutta  algorithm.  Complete  simulations  are  very  fast  on  a 
modern  computer. 

Equations,  parameters  and  initial  conditions  presented 
in  this  chapter  are  for  the  refined  version  of  the  CRM,  which  has 
evolved,  in  part,  from  the  studies  discussed  herein. 


* Effects  of  this  model  limitation  on  rise  simulation  of  very 
large  clouds  are  discussed  in  Chapter  4. 


2.2  SIMOLIFIED  EQUATIONS 

The  CRM  conservation  equations  are  those  developed  by 

(1  2) 

Huebsr.h  for  a water-surface  burst'  * ' as  modified  by  Norment  and 
Woolf  (ref.  3,  Appendix  C.l).  The  simplified  equations  given  here 
are  for  a pure-air  cloud  In  a dry  atmosphere.  The  more  complex 
equations  that  account  for  the  presence  of  water  and  soil  are 
given  in  ref.  3.  For  all  CRM  calculations.  Including  those  described 
here,  the  complete  equations  are  used.  Effects  of  water  and  soil  are 
deleted  solely  for  this  presentation  to  reduce  algebraic  complexity. 
To  further  increase  understanding,  each  equation  term  is  Identified 
as  to  function.  (Symbols  are  defined  in  Appendix  A below.) 

2.2.1  Momentum 


Terms  (a),  (b)  and  (c)  represent  forces  due  to  buoyancy,  eddy-viscous 
drag  and  entrainment  drag  respectively.  This  Is  the  only  equation 
that  differs  from  the  most  recent  predecessor  model  described  In 
ref.  3.  The  asymmetric  entrainment  factor  (1  - u)  and  the  factor  con- 
taining the  virtual  mass  have  been  dropped.  Deletion  of  these 

factors  causes  negligible  change  in  prediction  results. 


2.2.2  Center  Height. 


dz  _ 

at  ‘ 


u 


(2) 


2. 2. 3  Temperature 


dT 

3t  ■ 


(a) 


(b) 


(c) 


dm 

at 


(3) 
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Terms  (a),  (b)  and  (c)  account  for  the  effects  of  adalbatic  expansion, 
entrainment  and  dissipation  o-^  turbulent  eneryy  to  heat  (see  eq.  (4a)) 
respectively. 


2.2.4  Turbulent  Kinetic  Energy  Density 


dE 

dt 


(a)  (b) 

2 2 

T u V u dm 

T H.  m 3Y 
e c 


(c)  (d) 

1 dm 

m dt  ~ ^ 


(4) 


Terms  (a)  and  (b)  account  for  turbulent  energy  generated  by  the  mean 
flow  via  eddy-viscous  drag  and  momentum-conserving  inelastic-collision 
entrainment  i^espectively.  Term  (c)  represents  entrainment  dilution ^ 
and  (d)  represents  dissipation  to  heat,  where  the  turbulent  dissipation 
rate,  e,  is 


k,  (2E) 


3/2 


e = 


2.2.5  Mass 


dm  „ T 
dt  " T. 


(a) 

(b) 

S 

/ T 

Y uv  + ^ i 

(d) 


\ - 


(4a) 


(5) 


This  equation  is  based  on  observed  volumetric  growth  of  nuclear  clouds, 
and  term  (a)  contains  this  information.  The  other  terms  are  required 
to  convert  volumetric  to  mass  growth  riite.  Terms  (b)  and  (c)  account 
for  temperature  effects  caused  by  entrainment  and  conversion  of  tur- 
bulence energy  to  heat,  and  term  (d)  accounts  for  hydrostatic  expansion 
of  the  rising  cloud. 
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2.3  DIMENSIONLESS  PARAMETERS 


Three  empirical,  dimensionless  parameters  are  of  inter- 
est here.  All  are  held  constant  during  a cloud  rise  simulation. 

2. 3. 1 Energy  Conversion  Parameter,  k^. 

The  terms  in  eqs.  (1)  and  (4)  that  contain  kg  represent 
conversion  of  kinetic  energy  of  rise  to  turbulent  kinetic  energy. 
Since  I have  no  way  to  judge  a priori  what  the  value  of  kg  should 
be,  I have  taken  it  to  be  an  adjustable  parameter.  Good  prediction 
results  are  obtained  with 

kg  = max  [ 0.004,  min(0.1,  O.lW^^^)  | 


In  long  form,  this  equation  means: 


kg  = 0.1 

kg  = O.lW^''^ 
kg  = 0.004 

Restriction  to  a maximum  value  of  0.1  at  low  yields  was  found  to 
be  necessary  to  prevent  underprediction  for  low  yield  shots. 
Restriction  to  a minimum  value  of  0.004  at  very  high  yields  is 
discussed  in  sec.  4.2. 


W < 1 kT 

1 < W < 15,625  kT 
15,625  < W kT 


(6) 
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2.3.2  Turbulent  Energy  Dissipation  Parameter,  kj. 

kg  appears  in  eq.  {4a).  It  is  given  the  constant  value 
0.175.  which  is  unchanged  from  the  original  NRDL  version  of  the 
model 

2.3.3  Entrainment  Parameter,  u. 

The  entrainment  parameter  appears  in  eq.  (5).  As  ex- 
plained in  Appendix  C.l  of  ref.  3,  it  is  related  to  the  altitude 
gradient  of  linear  cloud  dimension  during  the  early  portion  of 
nuclear  cloud  rise.  In  the  preceding  CRM  version,  we  computed  m 
as  a function  of  yield  from  an  equation 

0 13 

M = 0.092W"-^'’ 

that  was  obtained  by  Norment  and  Woolf^^^  from  the  vertical  gra- 
dients of  horizontal  radii  for  a set  of  observed  doud  rise  data. 
Since  this  data  set  was  not  very  large,  and  since  one  cannot  be 
sure  that  horizontal  radius  is  the  proper  cloud  dimension  to  use, 

I have  simplified  this  relation  slightly.  Good  prediction  results 
are  given  by 

p = max[  max(0.12,  0.1W°*^)  , O.OIW^^^] 

In  long  form  this  equation  means: 

P = 0.12  , W ^ 6.192  kT 

p = 0.1W°-^  , 6.192  < W < 19,307  kT 

p = O.OIW^''^  , 19,307  < W kT 


(7) 
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It  was  found  necessci^^y  to  use  <i  -onstant  va'<ia  of  y - 0.12  for 
low  yield  shots  becajse  of  a tendency  to  cverpredict  for  low 
yields.*  The  change  from  1/10  to  1/3  yield  dependence  at  very 
high  yields  is  discussed  in  sec.  4.2. 


2.4  INITIAL  CONDITIONS 

Initial  time  (relative  to  detonation),  initial  gas 
phase  temperature  and  soil  loading  of  the  cloud  are  the  same  as 
in  the  original  DELFIC^^".  Only  initial  conditions  that  have 
been  changed  are  discussed  here. 

2. 4. 1 Explosion  Energy  Fraction  in  the  Cloud,  f. 

The  original  NRDL  water-surface  burst  model  used  an 
energy  fraction  of  33%.  For  the  first  DELFIC  simulations,  v^e 
increased  this  to  50%.  As  a result  of  extensive  simulations 
with  an  early  version  of  DELFIC,  we  have  used  f = 0.44W®*^^^ 

(W  in  kT),  which  is  such  a weak  function  of  yield  as  to  be  es- 
sentially constant. 

In  this  latest  version,  I simply  take  f = 0.45  for 

all  yields. 

2.4.2  I'emperature  of  Condensed- Phase  Matter. 

Fo*"  a surface  or  nesr  surface  burst,  a substantial 
amount  of  soil  which  does  not  reach  thermal  equilibrium  with 
the  fireball  gases  renains  in  the  solid  o*"  liquid  state.  The 
average  temperature  of  this  rrs.teria'(  3t  our  initial  time  is 


* If  the  yield  dependent  functions  for  k2  and  y are  allowed  to 
extend  into  the  low  yield  range,  the  tendencies  for  under  pre- 
diction because  of  high  k2  and  overprediction  because  of  low  y 
substentially  compensate.  However,  scatter  is  significantly 
greater  for  both  tops  and  bases  as  is  apparent  from  comparison 
of  observed  vs.  calculated  graphs  and  also  from  the  agreement 
statistics. 
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cor)j«»;tural  but  impertant.  since  the  energy  required  to  heet  it 
must  be  taken  from  that  available.  For  very  low  yield  shotSv 
the  soil  temperature  specified  by  the  original  DELFIC  Initial 
Conditions  tiodi'le  (Subroutine  TEMP)  has  been  found  to  be  too 
high^^^. 

I recommend  the  following  initial  soil  temperature 
for  future  use: 


Tj  ..  = 200  logjQ(W)  + 1000  , (8) 

where  W is  in  kT,  and  T.  . in  degrees  Kelvin. 

S 1 1 

2.4.3  Altitude^ 

In  the  preceding  CRM  version,  we  used  an  empirical 

equation  to  determine  cloud  center  height  relative  to  burst  point 

at  our  initial  time  (eq.  (1-20)  of  ref.  3).  I have  simplified 

1/3 

tiiat  expression  slightly  so  that  •'t  is  a function  of  W ' . Rood 
results  are  obtained  with 


^i  = 


“■GZ 


+ z 


HoB 


+ 90W 


1/3 


(9) 


2.4.4  Rise  Velocity. 

Since  1971  iry  version  of  the  CRM  has  used 

• (1" 

The  form  of  this  equation  resu!t,i  from  a simple,  approximate  anal- 
sis  of  initial  fireball  rise,  and  tSie  constant  2.2  is  chosen  to  fit 
observed  data. 
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2.4.5  Mass  and  Volume. 

The  energy  in  the  cloud  (45%  of  explosion  energy)  is 
used  to  heat  soil  and  air  to  their  respective  initial  tempera- 
tures. The  soil  inass  and  both  initial  temperatures  are  specified 
by  the  code^®^.  The  mass  of  air  is  computed  via  a simple  energy 
balance  as  described  in  ref.  3.  Cloud  volume  is  taken  to  be 
equal  to  the  volume  of  this  mass  of  air. 


2.5  CLOUD  SHAPE 

Observations  of  nuclear  clouds  show  that  vertical 
cloud  dimension  is  a simple  linear  function  of  height  during 
cloud  rise^^^.  This  holds  over  essentially  the  complete  rise 
history  of  nuclear  clouds  of  all  yields  for  which  we  have  data. 

In  the  preceding  CRM  version,  cloud  volume  was  computed  dynami- 
cally, and  vertical  cloud  radius  was  computed  as  a linear  function 
of  altitude.  The  slope  of  the  linear  function  was  taken  to  be 
y,  and  the  intercept  was  computed  for  an  ellipsoidal  cloud  of 
eccentricity  0.75  at  the  initial  time. 

Though  the  extensive  simulations  described  in  the 
next  chapter  have  shown  that  the  above  procedure  works  satis- 
fdctov'ily,  I have  deleted  it  from  this  version  of  the  CRM.  This 
was  done  to  avoid  the  possibility  that  this  rather  restrictive 
dependence  on  height  and  the  entrainment  parameter  y might  occa- 
S'^nally  produce  unrealistic  results. 

In  the  new  CRM,  the  clou:!  is  held  to  the  ellipsoidal 
shape,  eccentricity  0.75,  until  it  stops  rising,  after  which  its 
vertical  dimension  remains  fixed  while  volume  increase  is  accom- 
modated by  lateral  expansion.  The  eccentricity  value  0.75  is  an 
average,  with  standard  deviation  0.08,  found  by  Norment  and  Woolf 
for  ten  test  shots  that  cover  a yield  range  from  3.6  kT  to  15  MT^^\ 
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Norment  and  Woolf  found  little  variation  of  eccentricity  with 
height,  particularly  for  small  and  medium  yield  shots,  up  to  the 
altitude  at  which  stabilization  and  final  horizontal  expansion 
begins. 

For  a cloud  with  eccentricity  0.75,  H /R  = 0.66144. 


2.6  TERMINATION  OF  RISE  AND  EXPANSION 

Termination  switches  in  the  CRM  code  are  discussed  on 
pp.  48-51  of  ref.  3.  A normal  termination  was  via  the  R-RATE 
switch,  which  operates  when  the  relation 


0.014778 


1153 


(11) 


is  satisfied,  where  W Is  In  kT  and  At  In  seconds.  This  switch 
produces  reasonable  results  for  many  cases.  However,  it  fre- 
quently causes  trouble,  particularly  for  shots  at  the  extremities 
of  the  yield  range. 

For  low  yields,  the  shut-off  frequently  is  later  than 
desired  which  results  in  clouds  that  are  too  broad. 

At  the  other  extreme,  growth  rates  of  multi -megaton 
clouds  tend  to  oscillate  slowly,  and  If  a radial  expansion  rate 
trough  happens  to  dip  below  the  limit  given  by  eq.  (11),  shut- 
off occurs  prematurely. 

To  correct  the  high  yield  problem  I have  added  to 
eq.  (11)  the  requirement  that  the  cloud  rise  velocity  must  be  zero. 

For  low  yield  shots  I have  added  a new  switch,  the  U,EK 
switch.  This  terminates  cloud  growth  when  the  turbulent  kinetic 
energy  density  falls  below  a threshold.  Termination  occurs  when 
the  following  relations  are  satisfied 


E < max  1 10,  min(23  + 9 logj^;'. 


(12) 
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and 


u = 0 

In  long  form,  these  equations  mean  that  termination  occurs  when 

E < 10  . W f 0.0359  kT 

E < 23  + 9 logjQW  , 0.0359  5 W ^ 12,915  kT 

E < 60  , 12,915  < W kT 

and 

u = 0 

2 2 

Here  E is  in  units  of  Joules/kg  (i.e. , m /sec  ).  It  turns  out 
that  most  low  and  medium  yield  cases  are  terminated  by  this 
switch,  but  megaton  yield  cases  are  terminated  by  the  R-RATE  switch. 


3.  COMPARISON  OF  PREDICTED  WITH  OBSERVED  DATA 


3. 1 THE  OBSERVED  DATA  SET 

Almost  all  of  the  shots  that  I have  considered  in  this 
study  are  listed  in  Table  A. 2.1  of  ref.  3 (pp.  160-164).  For  each 
shot  the  data  consist  of:  yield,  ground  zero  height,  height  of 
burst,  observed  (at  or  near  shot  time)  altitude  profiles  of  pres- 
sure, temperature  and  relative  humidity,  observed  height  of  the 
stabilized  cloud  top  and,  in  most  cases,  observed  height  of  the 
stabilized  cloud  base. 

There  are  60  shots  in  the  set.  Yields  range  from  0.5  ton 
to  15  megatons.  In  recalibrating  and  validating  the  model  I have 
used  the  53  shots  with  yields  greater  than  0. D1  kT.  Observed 
heights  of  stabilized  cloud  bases  are  available  for  48  of  these. 

Top  and  base  heights  of  the  observed  stabilized  clouds 
are  plotted  against  yield  in  Fig.  1.  The  solid  lines  are  least 
squares  curve-  fits  to  the  data,  which  are  given  by 

(13) 

(14) 

where  Zt  and  Zo  are  in  meters  relative  to  burst  height  and  W is 
• * “ 
in  kilotons  . 

The  CRM  code  was  modified  such  that  it  could  batch  pro- 
cess the  complete  set  of  data,  print  individual  results,  print  sta- 
bilized height  agreement  statistics  and  plot  calculated  vs.  ob- 
served stabilized  cloud  heights.  Over  60  runs  were  made  as  I pro- 
gressed through  various  sensitivity  studies,  modifications  and 


Zj  = 3914W°*^^° 
Zo  = 1971W°*^°^ 


* These  equations  compare  well  with  their  counterparts,  eqs. 
(2.1)  - (2.6)  in  DASA-1251^^\  in  the  yield  range  0.1  to  100 
kT  (see  Table  1). 
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refinements  to  the  model.  The  complete  model.  Including  effects 
of  water  and  soil,  was  used  1n  all  calculations. 

3.2  COMPARISON  CRITERIA 

By  far  the  most  useful  criterion  for  judging  prediction 
adequacy  Is  comparison  of  calculated  with  observed  stabilized 
cloud  top  heights.  This  Is  true  for  the  following  reasons: 

• Stabilized  top  height  Is  the  most  critical  cloud 
property  In  determining  which  winds  are  Involved 
In  the  transport  and  deposition  of  fallout. 

• Stabilized  top  height  Is  by  far  the  most  accu- 
rately observed  cloud  property.  Other  stabilized 
cloud  properties,  such  as  base  height  and  radius, 
are  so  poorly  defined  as  to  be  virtually  useless 
In  assessing  simulation  accuracy. 

Therefore  I have  relied  mainly  on  comparison  of  stabilized  cloud  top 
heights  in  this  study. 

Of  course,  other  cloud  properties  are  of  consequence 
and  they  are  considered.  For  DELFIC  predictions  the  least  impor- 
tant is  the  height  of  the  base  of  the  stabilized  cloud  , since 
fallout  particle  cloud  bases  are  defined  by  the  altitudes  to  which 
the  particles  have  settled  at  stabilization  time.  As  regards  the 
stabilized  cloud  radius,  it  can  be  Important  in  determining  the 
crosswind  spread  of  the  fallout  pattern;  however,  assessment  of  the 
cloud  radius  prediction  accuracy  is  best  made  by  comparing  calculated 
with  observed  fallout  patterns,  which  is  not  considered  here.  Tem- 
poral development  of  nuclear  clouds  also  is  of  importance.  CRM  capa- 
bility to  simulate  all  of  these  cloud  properties  Is  analysed  in  the 
next  section. 

* This  corresponds  to  the  base  of  the  visually  observed  cloud  cap, 
which  is  also  the  quantity  plotted  In  the  b parts  of  the  figures 
herein. 
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3.3  RESULTS  FOR  THE  REFINED  MODEL 

3.3.1  Graphical  Comparisons  of  Stabilized  Cloud  Height. 

Calculated  vs.  observed  cloud  top  and  base  heights  In 
the  yield  range  0.01  - 15.000  kT  are  shown  In  Fig.  2.  The  tops 
comparison  Is  very  good  over  the  entire  yield  range.  For  the 
bases  there  Is  considerable  overprediction  for  low  yields,  whereas 
the  comparison  Is  good  for  medium  and  high  yields.  In  light  of 
the  scatter  of  observed  dcta  In  Fig.  1,  these  results  can  be  con- 
sidered to  be  excellent,  even  for  the  low  yield  bases. 

Calculated  vs.  observed  heights  for  a11  60  cases  (ex- 
tending the  yield  range  to  0.5  T)  are  shown  In  Fig.  3.  Obviously 
the  model  works  well  at  very  lo«f  yields  also. 


3.3.2  Stabilized  Cloud  Height  Comparison  Statistics. 

Agreement  statistics  used  In  the  refinement  arc  the  frac- 
tional RMS  deviation 


and  the  fractional  mean  deviation 
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Figure  2a.  (Calculated  vs.  observed  stabilized  cloud  top  heights 
(meters  relative  to  burst  height)  for  53  shots  In  the 
yield  rangt  0.021-15.000  kT  for  the  refined  CRM. 


23 


.00  4.00 


.00 


S.oo 


RUN-60 

3.00 


4.00 


LOGIO  OBS  CLOUD  BASE 


Figure  2b.  Calculated  vs.  observed  stabilized  cloud  base  heights 
(meters  relative  to  burst  height)  for  53  shots  In  the 
yield  range  0.021-15,000  kT  for  the  refined  CRM. 
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Figure  3a.  Calculated  vs.  observed  stabilized  cloud  top  heights 
(meters  relative  to  burst  height)  for  60  shots  In  the 
yield  range  0.0005-15,000  kT  for  the  refined  CRM. 
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Figure  3b.  Calculated  vs.  observed  stabilized  cloud  base  heights 
(meters  relative  to  burst  height)  for  50  shots  In  the 
yield  range  0.0012-15,000  kT  for  the  refined  CRM. 
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For  the  refined  model , we  have 


SHOTS  FOR  WHICH  W > 0.01  kT  COMPLETE  SET  OF  SHOTS 


TOPS 

BASES 

TOPS 

BASES 

FRMS 

0.15 

0.29 

0.14 

0.29 

FMD 

-0.009 

-0.16 

0.002 

-0.15 

Figure  1 and  eqs.  (13)  and  (14)  imply  that  much  of  tiie  scat- 
ter of  the  observed  data  can  be  accounted  for  by  expressing  Zy  and  Zg 
as  simple  power  functions  of  yield.  For  example,  the  least  squav'es 
(LS)  c.nalysis  by  which  eq.  (13)  was  determined  minimized  the  variance 

N , ,2 

,,  . EI'“9lO<^T>obs  - ’“9iO<^T>LsI  S % 

where  (Zy)L5  is  the  value  of  Zy  calculated  by  eq.  (13).  An  impc^tant 
criterion  of  the  value  of  the  CRM  is  to  compare  the  scatter  oT  obsfrrved 
data  around  the  least  squares  curve  with  the  scatter  of  observed  d?,ta 
around  the  CRM  calculated  results  Thus  the  result  of  eq.  (17)  is  to 
be  compared  with 


,,  - E('‘’9l0(^T>obs  - 'oSioU^jcRMf 

^ T ^CRM  N (18) 


where  (Zy)Qg|^  is  the  value  of  Zy  calculated  by  the  CRM. 
Results  of  this  comparison  are: 


TOPS 

0.0180 

0.00739 

53 

BASES 

0.0412 

0.0436 

48 
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For  the  tops,  the  CRM  affords  a proportional  reduction  in  variance 
of  59%  relative  to  the  simple  power  function  of  yield  (eq.  (13)). 

The  ratio  of  variances  ” 2.43)  can  be  used  to  test 

whether  or  not  the  reduction  in  variance  is  significant.  Using  the 
F distribution  test^  , the  reduction  is  significant  at  less  than  the 
2%  level.  (That  is,  the  probability  is  greater  than  98%  that  the 
reduction  in  variance  is  significant.)  For  the  cloud  bases,  on  the 
other  hsind,  the  power  function  of  y*eld  is  essentially  equal  to  the 
CRM  as  a predictor.  However,  as  noted  above,  cloud  base  heights  are 
very  poorly  defined  observational ly  and  are  not  critical  for  UELFIC 
fallout  predictions. 

It  is  important  to  realize  that  this  success  of  : ■ CRM 
is  by  no  means  its  only  virtue.  Of  more  importance  is  its  capability 
to  make  cloud  rise  predictions  under  atmospheric  conditions  that  are 
significantly  different  from  those  encounterod  at  the  Nevada  and 
Pacific  Test  Sites.  For  example,  tropopausa  heights  in  arctic  atmos- 
pheres or  in  northern  Europe  during  winter  are  much  lower  than  en- 
countered at  the  test  sites;  while  this  would  have  profound  effects 
on  local  distribution  of  fallout,  it  would  not  be  accounted  for  by 
use  of  the  empirical  power  function  equation  for  cloud  height. 

3.3.3  Cloud  Development  Comparisons. 

There  are  not  many  observed  cloud  rise  histories  avail- 
able that  are  detailed  enough  or  extend  to  early  enough  times  to 
be  of  value  here.  Moreover,  much  of  what  is  available  does  not 

correspond  with  our  data  set.  Figures  4 and  5 show  two  coinpari- 

* 

sons  available  from  the  data  on  hand.  Atmospheric  data  observed 


* It  is  instructive  to  note  that  the  observed  stabilized  cloud 
heights  shown  in  the  figures  do  not  correspond  closely  with  the 
values  used  in  our  comparisons  studies.  For  example,  the  observed 
U/K  Simon  top  height  in  Fig.  5 is  16  km,  whereas  the  value  used  in 
our  comparisons  is  13.3  km,  while  that  given  in  ENW^)  is  13.7  km. 
This  illustrates  the  considerable  variation  of  data  that  is  usual- 
ly found  in  different  sources. 
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ALTITUDE  (relative  to  MSL)  and  DIAMETER  (km) 


45  kT  tower 
;e  at  shot  t 


near  shot  time  are  used  in  the  simulations.  Figure  6 shows  cal- 
culated data  for  a 1.2  kT  surface  shot  with  a ground-zero  height 
of  1284  meters.  A standard  atmosphere  for  30°N,  January  was  used. 
Observed  stabilized  cloud  heights  for  Jangle-S  are  given  for  com- 
parison. 


3.3.4  Comparisons  of  Stabilized  Cloud  Dimensions  with  Standard  Equations. 

Table  1 lists  stabilized  cloud  data  for  sea-level  sur- 
face bursts  in  the  mid-latitude,  spring/fall  standard  atmosphere. 
Comparisons  with  results  calculated  by  eqs.  (13)  and  (14)  above, 
and  by  similar  equations  in  DASA-1251^^^  are  given,  including 
cloud  radii  comparisons.  Not  surprisingly,  agreement  of  CRM 
cloud  top  results  with  eq.  (13),  up  to  about  100  kT, 
is  usually  superior  to  agreement  with  the  DASA-1251  equations. 

On  the  other  hand,  for  higher  yields  the  reverse  is  true.  This 
latter  is  probably  a reflection  of  the  use  by  DASA-1251  of  a 
separate  equation  for  high  yield  shots  which  has  a lower  slope 
in  log-log  space  than  their  curve  for  low  yield  shots. 

For  low  yield  shots,  the  CRM  cloud  base  results  agree 
better  with  the  DASA-1251  equations  than  with  eq.  (14) , 

This  indicates  that  the  apparent  tendency  for  the  CRM  to  overpre- 
dict  cloud  base  heights,  and  produce  a cloud  that  is  too  thin, 
may  be  more  of  an  anomaly  of  the  data  set  used  here  than  a real 
probl em. 

The  radii  comparisons  in  Table  1 are  satisfactory  con- 
sidering the  element  of  subjectivity  that  is  inevitable  in  de- 
fining stabilization  time  both  in  the  field  and  in  simulations. 
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TABLE  1 

COMPARISON  OF  CRM  RESULTS  WITH  DASA-1251  EQUATIONS 


STABILIZED  CLOUD  DIMENSIONS  (meters) 


YIELD 

if 

TOP 

BASE 

JkTi 

SOURCE 

HEIGHT 

HEIGHT 

RADIUS 

10'^ 

CRM 

1010 

677 

265 

Eqs.(13)&(14) 

1130 

479 

DASA-1251 

1301 

722 

122 

10’^ 

CRM 

18C5 

1203 

460 

Eqs.(13)&(14) 

2103 

972 

DASA-1251 

2204 

1198 

327 

o 

o 

CRM 

3210 

2121 

844 

Eqs. (13)4(14) 

3914 

1971 

DASA-1251 

3734 

1987 

873 

10^ 

CRM 

6811 

4676 

1747 

Eqs. (13)4(14) 

7283 

3998 

DASA-125. 

6326 

3290 

2334 

10^ 

CRM 

12194 

7911 

4851 

Eqs. (13)4(14) 

13551 

8107 

DASA-1251 

14393 

9168 

6239 

10^ 

CRM 

18252 

10748 

14403 

Eqs. (13)4(14) 

25217 

16440 

DASA-1251 

21634 

13277 

16677 

10^ 

CRM 

32516 

16733 

39478 

Eqs. (13)4(14) 

46923 

33339 

DASA-1251 

32519 

19152 

44577 

10^ 

CRM 

59958 

17712 

178110 

Eqs. (13)4(14) 

87315 

67608 

DASA-1251 

48881 

27499 

119153 

* CRM  results  are  for  sea-level  surface  bursts  using  the  U.S, 
Standard  Atmosphere,  Mid-Latitude,  Spring/Fall.  DASA-1251 
height  results  are  computed  from  eqs. (2. l)-(2.6)  of  ref.  7, 
and  the  radii  are  computed  from  eq.(2.13}  of  ref.  7 for  a 
stabilization  time  of  ten  minutes. 
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3.3.5  Comparisons  Assessment. 

It  Is  clear  that  the  OELFIC  CRM  adequately  simulates 
the  general  features  of  nuclear  cloud  rise.  Moreover,  where 
observed  data  exist  to  provide  a basis  for  calibration  of  the 
model,  good  agreement  with  observation  is  obtained. 

3.4  SENSITIVITY  OF  RESULTS  TO  MODEL  STRUCTURE 

The  CRM  model  consists  of  twelve  coupled  differential 
equations  (eqs.  (1)  - (5)  plus  equations  to  account  for  water 
and  soil  content  of  the  cloud),  with  large  sets  of  Initial  con- 
ditions and  empirical  model  parameters.  The  model  Is  so  non- 
linear as  to  often  preclude  a priori  estlnatlon  of  the  effects 
of  changes  In  model  structure  or  Inputs,  unless  the  estimates 
are  based  on  extensive  experience  with  the  model.  In  this  sec- 
tion I present  some  additional  compa  Isons  results  and  discuss 
the  sensitivity  of  the  model  to  its  more  critical  or  controversial 
structural  aspects.  Sensitivity  to  Initial  conditions  and  atmos- 
pheric structure  Is  discussed  elsewhere^ 

3.4.1  DELFIC  CRM  circa  1970  and  1971:  Sensitivity  to  u^. 

Figure  7 shows  calculated  vs.  observed  stabilized  cloud 
tops  and  bases  for  the  version  of  the  CRM  described  in  ref.  3. 

The  major  fault  with  these  results  is  unexceptional  overpredict  lor. 
of  base  heights  for  low  yield  shots.  Otherwise,  the  results  are 
fairly  good  except  for  a tendency  to  underpredict  top  heights 
in  the  midyield-el uster  of  data. 

In  1971  the  initial  rise  velocity  was  changed  from  the 
power  function  of  yield  (eqs.  1.27  and  1.28  of  ref.  3)  to  that 
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Figure  7a.  Calculated  vs.  observed  stabilized  cloud  top  heights 
(meters  relative  to  burst  height)  for  53  shots  In  the 
yield  range  0.021-15,000  kT  for  the  ref.  3 version 
of  the  CRM. 
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Figure  7b.  Calculated  vs.  observed  stabilized  cloud  base  heights 
(meters  relative  to  burst  height)  for  48  shots  In 
the  yield  range  0.021-15,000  kT  for  the  ref.  3 
version  of  the  CRM. 
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given  by  eq.  (10)  above.  This  change  substantially  reduces  the 
Initial  rise  velocity  for  low  yield  shots. 

Figure  8 shows  results  for  the  CRM  that  are  Identical 
with  the  1970  version  except  for  this  change  In  Initial  rise 
velocity.  Comparison  of  Figs.  7 and  8 shows  that  reduction  of 
the  Initial  rise  velocity  results  In  an  Increased  tendency  to 
overpredict  both  tops  and  bases.  This  Is  a good  Illustration  of 
the  nonlinear  character  of  the  model.  It  also  Illustrates  rela- 
tively high  sensitivity  of  results  for  low  yield  shots  tc  changes 
In  Initial  rise  velocity. 

3.4,2  Sensitivity  to  p»  and  f. 

Prediction  results  are  critically  dependent  on  the 
values  of  the  turbulence  parameter,  k2.  and  the  entrainment  para- 
meter, M.  Increase  of  k2  causes  Increased  conversion  of  rise 
energy  to  turbulent  energy  (eq.  (4)),  and  Increased  turbulence 
drag  on  the  cloud  rise  (eq.  (1)):  these  effects  combine  to  pro- 
duce lower  cloud  heights.  Increase  In  y causes  Increased  entrain- 
ment of  air  (eq.  (5)),  which  In  turn  increases  entrainment  drag 
on  the  cloud  rise  (eq.  (1)),  reduces  cloud  temperature  (eq.  (3)), 
and  dilutes  turbulent  energy  density  (eq.  (4)).  The  overall  ef- 
fect is  that  lower  cloud  heights  result  from  increased  p. 

It  should  be  understood  that  low  yield  shots  are  more 
sensitive  to  such  changes  than  high  yield  shots.  The  reason  Is 
that  high  yield  clouds  are  brought  to  rest  In  the  stratosphere 
via  action  of  the  high  stability  encountered  there,  and  thus 
their  stabilized  heights  are  primarily  determined  by  the  height 
of  the  tropopause  and  structure  of  the  stratosphere. 

In  addition.  It  is  Important  to  know  that  effects  of 
changes  are  not  necessarily  similar  for  shots  of  similar  yields, 
and  large  changes  In  scatter  can  result  from  apparently  minor 
adjustments  of  parameters. 
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Figure  8a.  Calculated  vs.  observed  stabilized  cloud  top  heights 
(meters  relative  to  burst  height)  for  53  shots  In 
the  yield  range  0.021-15,000  kT  for  the  circa  1971  CRM. 
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Figure  8b. 


Calculated  vs.  observed  stabilized  cloud  base  heights 
(meters  relative  to  burst  height)  for  48  shots  in  the 
yield  range  0.021-15,000  kT  for  the  circa  1971  CRK. 
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As  shown  by  eqs.  (6)  and  (7}i  k2  and  y are  functions  of 
yield  over  most  of  the  yield  range  of  interest.  Extensive  calcu- 
lations, as  illustrated  by  Fig.  9,  have  shown  that  it  is  not  pos- 
sible to  adequately  match  calculated  with  observed  resul  if 
single,  yield  independent  values  are  used  for  these  parameters. 

For  airbursts,  stabilized  cloud  height  is  remarkably 
insensitive  to  change  in  explosion  energy  fraction,  f,  in  the 
cloud.  Variations  in  f from  20  to  50^  produce  only  minor  effects 
on  cloud  top  and  base  heights.  For  surface  bursts,  on  the  other 
hand,  large  amounts  of  soil  must  be  heated  and  the  available 
energy  must  be  sufficient  for  this.  Thus,  the  interplay  of  ini- 
tial soil  temperature  and  f,  discussed  above  on  pp.  14,15  becomes 
important. 

3.4.3  Sensitivity  to  Cloud  Shape. 

At  each  time  step  differential  equations  are  integrated 
at  the  point  of  the  cloud  center,  with  atmospheric  properties 
being  taken  for  that  point.  Cloud  volume  is  computed  from  mass 
of  cloud  air  using  the  ideal  gas  law;  therefore,  the  choice  of 
cloud  shape  appears  at  first  to  be  somewhat  arbitrary. 

Sensitivity  studies  show,  however,  that  results  are 
surprisingly  sensitive  to  cloud  shape.  This  is  caused  by  the 
effects  of  changing  the  vertical  radius,  H , on  the  differential 
equations,  H appears  explicitly  in  eqs.  (1)  and  (4),  and  in 
eq.  (5)  the  factor  S/V  reduces  to  3/H^  for  an  ellipsoidal  cloud> 
Furthermore,  shape  interacts  with  parameters  y and  in  complex 
ways  such  that  large  variations  in  scatter  are  induced  by  minor 
variations  in  shape,  apparently  depending  on  the  values  of  these 
other  parameters. 

Thus,  one  must  proceed  with  caution  in  manipulating 
cloud  shape  since  the  consequences  are  much  more  than  simply 
geometrical. 
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Figure  9a.  Calculated  vs.  observet!  stabilized  cloud  top  heights 
(meters  relative  to  burst  height)  for  53  shots  In  the 
yield  range  0.021-15t000  kT  with  constant  k^  and  y. 
The  only  substantial  difference  betvfeen  this  run  and 
that  of  F1g.  2a  is  that  k2  and  y have  values  0.1  and 
0.12  for  all  yields. 
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Figure  9b.  Calculated  vs.  observed  stabilized  cloud  base  heights 
(meters  relative  to  hurst  height)  for  48  shots  In  the 
yield  range  0.021-15,090  kT,  The  only  substantial 
difference  between  this  run  and  that  of  Fig.  2b  is  that 
ko  and  u have  values  0.1  and  0.12  for  all  yields. 


42 


g?.00  3.00  4.00  S.gp 


4.  EXTENSION  TO  VERY  HIGH  YIELDS 


4.1  LIMITATIONS  OF  THEORY  AND  DATA 

As  already  noted*  the  CRM  differential  equations  are 
integrctad  at  each  time  steo  at  the  point  of  the  cloud  center, 
using  atmospheric  conditions  at  the  point.  Cloud  volume  is  com- 
puted from  its  air  mass  using  the  ideal  gas  law.  Properties 
computed  at  the  center  are  taken  to  be  uniform  throughout  the 
cloud  volume. 

Even  for  low  yield  cases,  the  vertical  span  of  a cloud 
is  sufficient  to  cover  a substantial  range  of  atmospheric  condi- 
tions. For  large  yield  shots  this  span  of  conditions  is  immense. 
For  example,  at  our  initial  time  a 30  MT  surface  burst  cloud 
with  center  at  2.8  km  is  computed  to  have  a vertical  diameter  of 
4.4km.  At  stabilization  (using  the  15°N  annual  standard  atmos- 
phere) its  center  is  computed  at  36  km  with  vertical  diameter 
26  km.  For  a 100  MT  shot,  initially  the  cloud  center  is  com- 
puted at  4.2  km  with  vertical  diameter  6.6  km,  while  at  stabili- 
zation its  center  is  at  39  km  with  vertical  diameter  42  km. 

The  depth  of  the  troposphere  varies  from  roughly  10  - 
15  km.  Thus,  at  our  initial  time,  the  clouds  may  span  from  h 
to  mors  than  Hi  of  the  entire  troposphere.  The  100  MT  stabilizeo 
cloud  is  computed  to  be  distributed  through  the  entire  depth  of 
the  stratosphere,  and  even  overlap  into  the  mesosphere.  Repre- 
sentation of  such  huge  volumes  by  a point  in  the  calculations 
stretches  the  credibility  of  the  results. 

On  close  examination  of  results  for  very  high  yield 
shots,  one  does  indeed  find  non-physical  behavior  that  is  attrib- 
utable to  this  model  deficiency.  For  example,  extremely 
low  cloud  temperatures  (T  < 100°K)  are  computed,  and  the  cloud 
may  actually  lose  mass  for  a short  interval  of  time.  This 
occurs  because  the  factor  S/V  in  eq.  (5),  which  is  proportional 


to  the  ratio  of  cloud  surface  area  to  volume,  necessarily  becomes 
quite  rmall  for  very  large  clouds.  This  results  in  decreased  en- 
trainment, which  then  increases  velocity  via  reduced  entrainment 
drag  in  eq.  (1),  which  in  turn  acts  via  eq.  (3)  to  increase  tem- 
perature decay,  etc.  Obvious  abnormal  results  are  obtained  and 
sometimes  the  calculation  "crashes". 

Of  course,  these  problems  occur  to  lesser  extent  for 
smaller  clouds.  However,  for  smaller  yields  (W  - 15  MT)  we  can 
use  observed  data  for  calibration  and  validation.  It  is  important 
to  realize  that  the  capability  of  the  CRM  depends  as  much  on  this 
validation  as  on  the  theoretical  bases  of  the  model.  Consequent- 
ly, in  the  yield  range  from  20  - 100  MT,  prediction  results 
should  be  considered  as  conjectural.  The  user  must  realize  the 
theoretical  shortcomings  of  the  model  when  applied  to  very  high 
yields,  and  he  must  realize  that  above  15  MT  there  has  been  no 
validation  of  results. 

More  satisfactory  numerical  simulations  of  large  nucle- 
ar clouds  would  be  provided  by  two-  or  three-dimensional,  time- 
dependent  hydrocodes,  such  as  the  SHELL  code  used  by  the  Air 
Force  Weapons  Laboratory.  However  such  codes  have  their  own  dis- 
advantages, not  least  of  which  is  that  they  require  large  amounts 
of  computing  time  on  even  the  most  advanced  computers.  The  one- 
dimensional code  used  in  the  CRM  is  ideal  for  fallout  prediction 
needs,  except  as  noted  for  extremely  high  yields. 

4.2  ADJUSTMENT  OF  MODEL  PARAMETERS 

The  yield  dependent  functions  for  k2  and  u given  by  eqs. 
(6)  and  (7)  were  determined  such  as  to  produce  best  agreement  be- 
tween calculated  and  observed  stabilized  cloud  top  and  base 
heights  for  the  53  observed  cases  described  in  sec.  3.1,  which 
extend  to  15  MT  yield.  When  simulations  first  were  attempted  in 
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the  yield  range  from  20  - 100  MT,  I encountered  a tendency  for  the 
cloud  to  "run  away",  and  ultimately  for  the  calculation  to  "crash". 
The  main  cause  of  this  is  as  explained  in  the  preceding  section.  In 
addition,  the  low  values  assigned  k2  for  high  yields  decrease  eddy- 
viscous  drag  in  eq.  (1),  which  compounds  the  problem.  A straight- 
forward remedy  is  to  increase  both  and  u for  very  high  yield 
cases.  Thus,  a lower  limit  of  0.004  was  given  to  kg,  and  above  19  MT 
the  iield  dependence  of  u was  changed  from  to  These 

adjustments,  which  are  expedient  rather  than  empirical,  enhance  the 
probability  that  simulations  will  run  to  completion  for  yields  up 
to  100  MT.* 


* There  can  be  no  guarantee  that  "plausible"  results  will  be  ob- 
tained. For  example,  apparently  minor  defects  in  atmospheric 
data  can  cause  trouble.  Relative  humidity  data  usually  are  not 
given  at  altitudes  above  1C  km  in  standard  atmosphere  tabula- 
' ‘ 'ns.  relative  humidity  is  taken  to  be  zero  above  10  km, 
prema";  r'fc  termination  of  cloud  rise  may  occur  in  the  yield 
range  50-100MT.  This  has  been  corrected  by  simply  assuming 
reasonable  values  of  relative  humidity  above  10  km. 
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APPENDIX  A 
GLOSSARY  OF  SYMBOLS 

(mks  units  are  used  except  where  noted  otherwise) 

specific  heat  of  air  at  constant  pressure 

differential  operator 

cloud  turbulent  kinetic  energy  density 

fraction  of  the  total  explosion  energy  in  the 
cloud  at  the  initial  time 

gravity  acceleration  constant 
vertical  cloud  radius 

rise  kinetic  energy  to  turbulent  energy  con- 
version parameter  (dimensionless) 

turbulent  energy  dissipation  constant  (di- 
mensionless) 

cloud  mass 

cloud  virtual  mass  at  the  initial  time 
ideal  gas  law  constant  for  air 
horizontal  cloud  radius 

time 

cloud  temperature  (°K) 
ambient  temperature 

temperature  of  soil  matter  in  the  cloud 


\ 
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U cloud  center  velocity 

y mBx(u»  \/5F) 

1' 

V cloud  volume 

W explosion  energy  yield  (kT) 

z vertical  coordinate 

G turbulent  kinetic  energy  density  dissipation 

rate  In  the  cloud 

u entrainment  parameter  (dimensionless) 
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APPENDIX  B 


REVIEW  OF  THE  EUCLID  RESEARCH  GROUP  ANALYSIS 
OF  THE  ^LOUD  RISE  MODULE 


B.l  INTRODUCTION 

In  an  effort  to  remedy  certain  problems  they  were  having 
with  the  DELFIC  CRM,  the  Ballistics  Research  Laboratories  retained 
Euclid  Research  Group  (ERG)  to  analyse  the  CRM  model  described  in 
ref.  3,  and  to  recommend  solutions  to  the  problems.  A report 
that  contains  the  ERG  analysis  and  recommendations  has  recently 
been  published^^^^. 

Those  problems  which  in  fact  do  exist  are  discussed  and 
explained  in  the  body  of  this  report.  The  purpose  of  this  review 
is  to  show  the  consequences  of  following  the  ERG  recommendations, 
and  to  point  out  why  these  recommendations  fail. 

B.2  THE  ERG  ANALYSIS 

The  ERG  analysis  consists  of  three  major  parts: 

• Consideration  of  the  effects  on  energy  transfer 
balance  of  two  extraneous  factors  in  the  momen- 
tum equation  (eq.(l)). 

• An  analysis  of  the  entrainment  equation  (eq.(5)). 

• Discussion  of  cloud  shape. 

Each  of  these  is  discussed  below. 
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B.2.1.  The  Momentum  Equation . 


The  momentum  equation  given  In  ref.  3 (eq.(l))  has  two 

extraneous  factors,  — 7 — • which  ERG  labels  the  "virtual  mass 

m + m^ 

factor"  (VNF),  and  (l-u)  which  ERG  Inappropriately  labels  the 
"shape  factor"  (SF).  The  VMF  factor  Is  a relic  of  the  original 
version  of  the  CRM  as  developed  by  the  Naval  Radiological  Defense 
Laboratory  (NRDL).  The  SF  factor  was  Introduced  by  n^self  In  the 
ref.  3 version  of  the  CRM. 

By  means  of  a theoretical  analysis,  ERG  shows  that  these 
factors  cause  imbalances  In  energy  transfer  In  the  cloud.  Fortu- 
nately, It  turns  out  that  the  factors  have  negligible  effect  on  the 
cloud  rise.  They  have  been  deleted  from  the  refined  CRM. 


B. 2. 2 The  Entrainment  Equation. 

A key  feature  of  the  NRDL  cloud  rise  model,  developed  by 
Huebsch^^^,  was  the  formulation  of  the  entrainment  equation.  The 
success  of  the  model  was  to  considerable  extent  due  to  this  formu- 
lation. In  Its  original  form,  the  equation  was 


Some  years  ago,  I derived  this  form  of 
the  equation  by  combining  the  equation  for  mass  conservation  with 


* See  the  last  paragraphs  In  secs.  3.2.3.  and  4.2.3.  of  the  ERCi 
report. 
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where  X was  a constant. 


a simple  function  for  vertical  gradient  of  linear  cloud  di- 
mension that  was  obtained  from  observed  data.  This  analysis 
showed  that  additional  terms  are  required  In  the  entrainment 
equation  (see  Appendix  C.l  of  ref.  3),  and  the  CRM  model  was 
revised  accordingly. 

ERG  devotes  considerable  effort  to  rationalizations  of 
why  these  changes  In  the  entrainment  equation  cause  problems 
with  the  simulations.  The  problems  discussed  In  this  connection  are 
those  that  affect  extremely  high  yield  cloud  simulations. 

These  problems  arise  from  the  mathematical  representation  of 
the  entire  cloud  volume  by  a point,  a fundamental  limitation 
of  the  model  that  Is  discussed  In  Chapter  4 above. 

In  fact  the  revised  entrainment  equation  Is  theoreti- 
cally sound.  This  Is  substantiated  by  the  ERG  energy  trans- 
fer analysis,  which  shows  tot.al  energy  balance  when  the  VMF 
and  SF  factors  are  removed.  Furthermore,  the  form  of  the 
equation  Is  solidly  based  on  observed  data.  And  finally,  the 
appropriateness  of  the  revised  formulation  is  supported  by 
the  success  of  the  model  in  sin%tlat1ng  observed  cloud  behavior, 
as  Is  well  documented  in  the  body  of  this  report. 
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B.2.3  Cloud  Shape. 


As  implied  in  sec.  2.5  above,  the  formulation  of  cloud 
shape  in  the  context  of  this  model  has  always  been  a minor 
problem.  The  formulation  used  in  the  ref.  3 version  of  the 
CRM  is  discussed  in  sec.  4.2.4  of  the  ERG  report.  I concur 
that  that  approach  is  less  flexible  than  is  desirable,  and 
the  possibility  exists  that  unrealistic  results  may  occur, 
though  it  does  not  happen  in  simulations  of  test  shots.  The 
new  CRM  uses  an  ellipsoidal  cloud  with  eccentricity  0.75  as 
recommended  by  ERG.  However  this  choice  is  not  made  for 
reasons  of  geometry  alone,  which  is  the  sole  basis  of  the 
ERG  discussion. 

It  turns  out  (see  sec.  3.4.3  above)  that  the  overrid- 
ing effect  of  manipulating  cloud  shape  is  to  radically  alter 
the  rise  via  the  effect  of  the  vertical  cloud  radius,  H^,  on 
the  dynamics.  Thus,  the  choice  of  an  ellipsoidal  cloud  was 
made  on  the  basis  of  comprehensive  simulation  results,  and 
the  particular  eccentricity  value  was  chosen  to  conform  to 
observed  cloud  behavior^^^. 

B.3  THE  ERG  RECOMMENDATIONS 

ERG  makes  the  following  recommendations  for  immediate 
implementation: 

• Return  the  momentum  and  entrainment  equations  to 
the  form  used  in  the  original  NRDL  version. 

• Use  an  ellipsoidal  cloud,  eccentricity  0.75,  un- 
til the  cloud  passes  the  tropopause  or  stops 
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rising,  after  which  volume  growth  is  accommo- 
dated by  lateral  expansion.  (The  NRDL  version 
used  this  same  algorithm  except  that  a spheri- 
cal , instead  of  ellipsoidal,  cloud  was  used.) 

• Set  the  parameters  k2  and  u to  values  of  0.1 
and  0.25.  (The  NRDL  version  used  values  of 
0.1  and  0.2  for  these  parameters.) 

Implementation  of  these  recommendations  yields  the 
results  shown  in  Fig.  B.l  — an  obvious  disaster.  (See  sec.  3.1 
above  for  a description  cT  the  observed  data  set.)  This  may 
seem  surprising  since  the  model  is  so  close  to  the  one  developed 
by  Huebsch  when  he  worked  at  NRDL.  However,  NRDL  never  performed 
a satisfactory  validation  study  of  their  model,  apparently  being 
satisfied  with  results  that  were  merely  plausible. 

B.4  CONCLUSIONS 

The  ERG  recommendations  amount  to  a return  of  the  model 
substantially  to  its  form  of  twelve  years  ago,  as  originally  devel- 
oped at  the  Naval  Radiological  Defense  Laboratory  for  water  surface 
bursts.  As  is  shown  by  Fig.  B.l,  that  version  of  the  model  is 
still  unacceptable. 

The  ERG  analysis  is  purely  theoretical;  no  sensitivity 
computations  were  made  to  check  conclusions  arrived  at  by  compli- 
cated rationalizations-  Likewise,  no  attempt  was  made  to  deter- 
mine the  effects  of  recommended  changes  on  model  performance  over 
the  complete  range  of  yields  of  interest.  In  light  of  the  non- 
linearity of  results  obtained  by  numerical  solution  cf  the  twelve 
coupled  differential  equations  that  form  the  basis  of  the  modal, 
there  is  small  wonder  that  the  analysis  failed  in  its  objective. 
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Figure  6.1a.  Calculated  vs.  observed  stabilized  cloud  top  heights 
(meters  relative  to  burst  height)  for  53  shots  In 
the  yield  range  0.021-15,000  kT  for  the  CRM  version 
recommended  by  Euclid  Research  Group. 
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Figure  B.2b.  Calculated  vs.  observed  stabilized  cloud  base  heights 
(meters  relative  to  burst  height)  for  48  shots  In  the 
yield  range  0.021-15,000  kT  for  the  CRM  version  re- 
commended by  Euclid  Research  Group. 
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